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Abstract 

We present a quasiclassical theory of a decay accompanied by bremsstrahlung with a special 
emphasis on the case of ^^'^Po, with the aim of finding a unified description that incorporates 
both the radiation during the tunneling through the Coulomb wall and the finite energy of the 
radiated photon up to ~ Qa/y^, where Qa is the Q-decay Q-value and r] is the Sommerfeld 
parameter. The corrections with respect to previous quasiclassical investigations are found to 
be substantial, and excellent agreement with a full quantum mechanical treatment is achieved. 
Furthermore, we find that a dipole-quadrupole interference significantly changes the a-7 angular 
correlation. We obtain good agreement between our theoretical predictions and experimental 
results. 
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I. INTRODUCTION 

A characteristic feature of the a decay process is the quantum mechanical tunnehng |l| 
through the so-called Coulomb wall generated by the electrostatic interaction of the a par- 
ticle with the constituent protons of the daughter nucleus. Bremsstrahlung in a decay is 
intriguing because of the classically incomprehensible character of radiation emission during 
the tunneling process. Considerable attention has therefore been devoted to both experi- 



tne tunneiingprocess. uonsideraDle attention nas tnereiore been devc 
mental [2I, U Ij, Q] as well as theoretical investigations 0, I^, I^, 1^, 
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with the aim of elucidating the role of tunneling during the emission process. It is neces- 
sary to emphasize, however, that the term "radiation during the tunneling process" has a 
restricted meaning as the wavelength of the photon is much larger than the width of the 
tunneling region and even larger than the main classical acceleration region. It is therefore 
not possible to identify the region where the photon was emitted. Besides, it is possible to 
write the matrix element of bremsstrahlung in different forms using operator identities. As 
a result, the integrands for the matrix element, as well as the relative contributions of the 
regions of integration, will be different depending on the operator identities used, although 
the total answer remains, of course, the same. This was demonstrated, e.g., by Tkalya in 
Refs. [ll,[3. 

In the present paper, we revisit the theory of bremsstrahlung in the a decay of a nucleus 
with a special emphasis on the quasiclassical approximation. The applicability of this ap- 
proximation is ensured by the large value of the Sommerfeld parameter rj (see below). We 
investigate the range of validity of the result obtained by Dyakonov and show that it is 
restricted by the condition x ^ 1/y^ where x = E^/Qa (here, is the photon energy, and 
Qa is the a-decay Q- value). Our quasiclassical result has no such a restriction although we 
assume x <^ 1. It is consistent with the results of both Dyakonov and Papenbrock and 
Bertsch js] in limiting cases. For the experimentally interesting case of the a decay of ^^''Po, 
our result is valid with high accuracy up to x ~ 0.1. 

Another subject investigated here is the angular distribution of emitted photons. The 
a particle, initially in an 5* state, may undergo a dipole transition to a P final state, or a 
quadrupole transition to a D state. While the quadrupole contribution is parametrically 
suppressed for small photon energies, the effective charge prefactor for the quadrupole con- 
tribution is large. The dipole-quadrupole interference term vanishes after angular averaging. 
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but gives a significant contribution to the differential photon emission probabihty, resulting 
in a substantial deviation from the usually assumed dipole emission characteristics. 

Very recently, the results of our high-statistics measurement of bremsstrahlung emitted 



16|. Due to the limited solid-angle 



in the a decay of ^^°Po have been published, see Ref. 
coverage of the detectors used in this experiment, it was necessary to account for the a-7 
angular correlation. Taking into account the contributions of the dipole and quadrupole 
amplitudes in the data analysis, as derived in the present paper, good overall agreement 
between theory and experiment is observed. 

This paper is organized in four sections. In Sec. [Tll we investigate the leading dipole 
contribution to the differential bremsstrahlung probability and evaluate the corresponding 
amplitude in the quasiclassical approximation. The quadrupole contribution to the ampli- 
tude and its interference with the dipole part is analyzed in Sec. IIIII Conclusions are drawn 
in Sec. HVl Two appendices provide details on the methods used in the calculations. 



II. DIPOLE EMISSION 
A. Emission Probability 

It was shown in Ref. 8|] that the differential bremsstrahlung probability dP/dE^ as a 
function of the energy of the radiated photon in the dipole approximation has the form 



dP 4e2z2 



cff I A 



where natural units with h = c = = 1 are applied throughout the paper, e is the proton 
charge and /i is the reduced mass of the combined system of a particle and daughter nucleus, 
V = V{r) = ziZ — z)a/r is the potential of the daughter nucleus felt by the a particle. The 
functions Ri and Rf are the radial wave functions of the initial and final states corresponding 
to the angular momenta / = and / = 1, respectively (see App. [B]). The effective charge 
for a dipole interaction between an a particle with charge number z = 2 and mass number 
4 emitted from a parent nucleus with charge number Z and mass number A is (see also 
App.E]) 

^cff - ^cff - ^ ' l^j 

where the latter value is relevant for the experimentally interesting case of the a decay of 
^^°Po [Z = 84). Evaluating the effective charge with accurate values for the masses of the 



alpha particle and the daughter nucleus (^"^Pb, Z = 82), as given in Ref. 17| yields a value 
of Z^g^ = 0.399. 

In the present paper, we calculate the matrix element Ai in the quasiclassical approxi- 
mation taking into account the first correction of the order 0{ri~^), and the corresponding 
result is given below in Eq. (jS]). However, before presenting and discussing our formula for 
the matrix element Ai, let us briefiy review several results for Ai obtained earlier in the qua- 
siclassical approximation. These are illustrative with respect to their range of applicability 
and with respect to the importance of the tunneling contribution. 



B. Dipole Transition Matrix Element 

Various approximations have been applied for the evaluation of the matrix element Ai 



in Eq. ([T]) [8|, |9|, [lOl, [ill, ll^. The approximations are intertwined with the identification of 
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particular contributions to the real and imaginary parts of the matrix element At due to 
"tunneling" and due to "classical motion" of the a particle. 

We use the convention that the complex phase of the matrix element Ai should be chosen 
in such a way that it becomes purely real in the classical limit -E^ — > 0. Our definition of 
Ai is consistent with that used in Ref. and differs by a factor i from the definition used 



a 



in Ref. 

Equation (5) in the work of Papenbrock and Bertsch 8[ contains a fully quantum mechan- 
ical result for the photon emission amplitude Ai, expressed in terms of regular and irregular 
Coulomb functions, without any quasiclassical approximations. However, the physical in- 
terpretation of this result depends on a comparison with a quasiclassical approximation, 
as only such a comparison clearly displays the importance of the finite photon energy and 
the emission amplitude during tunneling. Papenbrock and Bertsch js] therefore present and 
discuss a quasiclassical expression for the imaginary part of their matrix element (real part 
for our convention), ignoring the contribution from the tunneling process to the emission 
amplitude. Note that the quasiclassical expression of Papenbrock and Bertsch provides 
a very good approximation for the imaginary part of their matrix element up to very large 
photon energies with x ~ 0.6. 

n 

In contrast, the quasiclassical result of Dyakonov [9j] is valid only for very small photon 
energies (x <^ 1/y^), but includes contributions from tunneling. Here, we unify the treat- 
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ments of Refs. j^, and obtain a quasiclassical differential emission probability dP/dE^, 
which includes the effect of photon emission during the tunneling process and which is 
substantially more accurate for higher photon energies than that of Dyakonov 

The quasiclassical approximation for the wave functions of the system of an a particle and 
a daughter nucleus in the initial and final states is valid for large values of the Sommerfeld 
parameters rjij = z{Z — z)e^^/kij with ki = \/2^Qa and kf = ^j2^{Qa — Ej). The indices 
i and / are reserved for initial and final configurations throughout the paper. The value of 
rji for the decay of ^^'^Po, which is the experimentally most interesting nucleus, is 22.0, while 



the Q-value is Qa = 5.40746 MeV 



17|. 



If one neglects the contribution to the matrix element Ai from the region r < Tq, where Tq 
is the nuclear radius, then one is consistent with the simple assumption for the potential as a 
square well for r < tq and a pure repulsive Coulomb potential for r > tq. Implementing this 
procedure according to Papenbrock and Bertsch , one obtains the following approximation 
for the real part M of the matrix element of Ai, 



M = ReM 




2 ki 

71 kf 

°° dr 

— Fi{r]f,kfr) Fo{r]i,kir) . (3) 

Here, Fq and Fi are the regular Coulomb radial wave functions corresponding to angular 
momenta / = and / = 1, respectively. The importance of the contribution of the region 
^ ^0 was discussed in Refs. jsl,]^!. For relatively small photon energies, which are interesting 
from an experimental point of view, this contribution is not significant, and we do not 

or high photon energies, the contribution of the region 



13|). The explicit form of M is given by Eqs. (6) and (7) 



consider it in the present paper, 
r ~ rg can be important (see Ref. 
of Ref. BQ. 

We have calculated the integral in Eq. using quasiclassical wave functions, keeping the 
correction of order 0{f]~^) and ignoring terms of order and higher in the expansion 

for large "mean" Sommerfeld parameter f/ = {r]i + rif)/2. The result of such a quasiclassical 
calculation for the real part of Ai reads 



M 



X ^e 



2 ki ki kf Tji 
71 kf ki + kf fj 

' ' (4) 
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FIG. 1: The ratio of M = ReTW, see Eq. ([3]), to various approximations, as a function of x = E^/Qa 
for for the case of ^lOpo (r/j = 22.0). The soHd hue shows M/M with M taken from Eq. ([3D, 
illustrating the excellent agreement of the quasiclassical matrix element with the exact result in 
the range x < 0.3. For the dashed line we use , which is obtained from M by omitting the 
correction of order ©(ry"^), showing a deviation of less than 5%. For the dash-dotted line the 
asymptotics M"^ as given in Eq. ([5]) is used. The deviation of the latter curve from the others 
illustrates that the "x 0" -asymptotics is indeed only valid for x <^ 



where ^ = rjf — rji, and Ka{h) is the modified Bessel function. The derivative is K'^{h) 



■^Kaib). Papenbrock and Bertsch |8j also calculated the real part of the integral ([3]) in the 
quasiclassical approximation. Note, however, that their term of order 0{f]~^) [see Eq. (14) 
of Ref. 81]] contains an additional factor -\/3/2 in comparison to Eq. (j4]). 

Figured] shows the ratio Re = M/ M as a function of x for rji = 22.0, corresponding 

to the case of ^^'^Po. Note that even at x = 0.3, the deviation of the quasiclassical result 
with the correction 0{f]~^) taken into account is less than 1% (solid line), while without 
this correction (M°), the deviation is about 5% (dashed line). 

Strictly speaking, the validity of the evaluation of the transition matrix element ([T]) with 
the wave functions taken in the quasiclassical approximation requires special consideration 
(see §51 of Ref. [3]) because of possibly noticeable contributions from the vicinities of the 
turning points. However, one can show that, for ^ <^ f] (or x <^ 1), the contributions of 
the vicinities of the classical turning points = 2rii/ki and = 2rif/kf in our case) are 
small. For x < 1, these contributions are no longer negligible. 

For X <^ 1, we have ^ = rji (x/2 + 3x^/8 + ...). If x^?7j -C 1 (even if xrji < 1), we can 
replace in (j4]) ^ by xrii/2 and make the substitution f]f ^ f]i and kf ^ ki. As a result, we 
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obtain the following asymptotics of M: 



2tx 
-Ki 



xrji 
2 



-iTxrii/4: 



ixm/2 



2 



ixTii/2 



2 



(5) 



From the dash-dotted line of Fig.[T], we see that the ratio M / M"^ deviates substantially from 
unity, illustrating that the applicability of Eq. ([5]) is indeed restricted to very small values 
of X. 
Since 

/■oo 

-ft'jiy(a^) = exp(— 7rz//2) / cos(x sinh t — z/t) c/t , (6) 

Jo 

the quantity without the Oijjl^) correction exactly coincides with the real part of the 
amplitude obtained by Dyakonov j^, where 



My) 



ki J f XTji 



X 



27r ^ V 2 
iyexp{-TTy) 

dt sinh(t) exp[iy (sinht — t)] 



(7a) 



(7b) 



showing that this result is applicable only to very small photon energies. Aij) contains both 
a real and an imaginary part and thus takes photon emission during tunneling into account. 
It was pointed out by Dyakonov that the imaginary part of A^d at xrji ~ 1 is of the same 
order as the real part. For instance, in the case of the a decay of ^^'^Po with rji = 22.0, we 
have x?7j ~ 1 for x ~ 0.05. This indicates that the imaginary part of TWd is important. 

Our quasiclassical approximation Ai for the dipole transition matrix element Ai from 
Eq. ([1]), which includes the contributions from the tunneling of the a particle (we would like 
to refer to this result as the "unified result" in the following sections of the current paper) 
has the form: 



M 



2ki kikf r], 



V 



n kf ki + kf 7] 
J{y) = i?/exp(-7ry) 

dt sinh(t) exp[iy {t — sinht)] , 

/•oo 

exp(— Try) / dt exTp[iy {t — sinht)] . 
Jo 



My) 



(8a) 

(8b) 
(8c) 
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FIG. 2: The ratio ImM/ReM for the case of ^^"Po {rji = 22.0). The graph illustrates that the 
imaginary part of the matrix element M is quite substantial even for moderate values of x. 

The derivation of Eqs. [8] (see App. [B] for more details) involves a shift of the integration 
region into the complex plane, which is needed to take the tunneling region into account in a 

nn 

quasi-classical treatment, as explained in [7|, |9|. Note that J{y) is the complex conjugation of 
the function J^iv) as defined in Eq. ([7]). Although this is irrelevant for the calculation of the 
bremsstrahlung emission probability, it is important for the dipole-quadrupole interference 
term discussed in Sec. IIIII The region of applicability of Eq. ([8]) is much wider than that of 
Eq. ([7]), because at small x there is no additional restriction x ^ 1/ which may otherwise 
constitute a strong limitation at large value of rj, as it was shown above. Moreover, Eq. ([8]) 
contains a correction 0{f]~^), which is also essential. 

Our unified result ([S]) can be used with high accuracy up to x ~ 0.3. The real part of Ai 
is identical to M given in Eq. (jl]) and is discussed already in detail above (see also Fig. [1]). 
Figure [2] shows the ratio ImAi/KeAi = ImAi/M as a function of x at ?7i = 22.0, i.e. for 
the case of ^^°Po. One can see that ImA^ is not small in comparison to ReA^. Thus, the 
imaginary part gives a noticeable contribution to dP/dE^ even for small x, and should not 
be neglected. This point was also emphasized in Refs. jol. [lol|. 



C. Quantitative comparison of various quasiclassical results 

As a last step, we compare in Fig. [3] the differential bremsstrahlung probability dP/dE^ 
for the bremsstrahlung accompanying a decay of ^^'^Po obtained with the use of our matrix 
element Ai given in Eq. the matrix element M [Eq. (jl])], corresponding to the quasi- 
classical approach of Ref. 8|], and the matrix element Aij) [Eq. ([7j)], to the full quantum 
mechanical formula given by Eq. (5) in Q]. A detailed comparison is shown in the bottom 
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FIG. 3: Differential brancliing ratio dP/dE^ for bremsstrahlung emission during a decay of ^Po 
in units of inverse keV (top panel) and relative to the fully quantum mechanical calculation (PB) of 
Papenbrock and Bertsch [8|] (bottom panel). The thick solid curve corresponds to our quasiclassical 
result, as given in Eq. ([8]). The other results are based on the quasiclassical treatment of Ref. Q] 
(PB), the semiclassical treatment of Dyakonov 9|], and on Eq. ([8]) of the present work but neglecting 
the correction terms of order 0{f]~^). 



panel of Fig. [3] in the 7 energy range < -E^ < 600 keV. While the result of Dyakonov 9| 
(dotted line) deviates roughly by a factor two at = 600 keV, our result (thick solid 
line) agrees with the exact quantum mechanical treatment within about 2 % at this photon 
energy; the inclusion of the 0{f]~^) correction is crucial in obtaining this agreement, as is 
evident from the supplementary curve in the bottom panel, where we omit the Ji term from 
Eq. (IScl) . The quasiclassical approximation of Papenbrock and Bertsch jsl, obtained by ne- 
glecting the imaginary part of the matrix element (dashed line), deviates by more than 15 % 
at = 600 keV. As expected, at low photon energies all results agree with each other. 
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FIG. 4: Interference term xi^)^ defined in Eq. (jlip . for the bremsstrahlung accompanying a decay 

of 210Po. 

III. QUADRUPOLE EMISSION 

We are now concerned with the quadrupole component of the bremsstrahlung probabihty 
and the angular distribution of the radiation due to interference with the dipole components. 
The outgoing a particle defines an axis of symmetry. We therefore we may use dQ = 
27T sin 6d6 in order to describe the solid angle element of the photon spanning an infinitesimal 
range of polar angles 9 with respect to the direction of the emitted a-particle. We assume 
that a summation with respect to photon polarization is performed. Within the dipole 
approximation, Eq. ([1]) gives rise to an angular distribution of the form 



dE^dVt 



7(1) 

"efr 



sin^^ 



dip ^r^l^E^ 



\M\ 



(9) 



where the index refers to the dipole approximation. Including the quadrupole term (see 
App. [B]), this formula should be generalized to 



d'^P _ sin^ 9 
d^P 



dE^dfl 



Z^^e'^^M cos 9 



dE^dVt 



+ X cos^] + C(Ar' 



(10) 



dip 



with 



X 



7(2) 



Re 



,i(5i-52) 



(11) 



where and 82 are the Coulomb phases corresponding to the angular momenta / = 1 and 
1 = 2, respectively. The effective quadrupole charge Z^^ is approximately given by (see 
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App.\E> 

Z'S . 2 + = 1.954 (12) 

for the case of ^^°Po, and is roughly five times larger than the dipole effective charge given 
in Eq. ([2]). Exact masses lead to the same result Z^^^ = 1.954 (up to the last decimal 
digit indicated). 

Calculating the quadrupole matrix element A/" within the quasiclassical approximation 
and keeping the leading in 1/f/ term, we obtain 



Af= - J^j^^vMO. (13a) 

We here neglect a parametrically suppressed next-to-leading order correction to the in- 
terference term, discussed in more detail in Appendix [B], and we introduce the nota- 
tion V = \/2Qal yU, where v approximately equals the final velocity of the a particle for 
bremsstrahlung emission with x <^ 1. For ^^''Po, we have v ^ 0.05. Note that for a large 
Sommerfeld parameter ?7j we have 

gi(52-5i) + A. (14) 



Because of this "i" in the right-hand side of Eq. f[Hj) . the imaginary parts of the functions 
J(.^) and Ji(0 become very important for the interference term. The quantity x(3^) defined 
in Eq. ( fTTj) . which determines the relative magnitude of the interference term, vanishes for 
a; — in the leading in 1/f/ approximation as 



TT Z^^^ 

X^i^^^-^ri.x. (15) 



For ^^°Po, this evaluates to x ~ 8.64 x. Therefore, the value of the coefficient xip^) becomes 
significant already at very small photon energies (x ~ 0.1 at E^^ ~ 0.06 MeV). The coefficient 
x[x\ calculated with our quasiclassical results M. and A/", is plotted in Fig. [Has a function 
of a; = E^/Qa for the experimentally interesting case of ^^°Po. 

When integrating Eq. (fTOl) over the total solid angle, the interference term drops out 
and the remaining quadrupole contribution to the differential emission probability is of the 
order 0{N"^)\ for the case of ^^°Po this contribution amounts to less than 1.5 % of the leading 
dipole term for photon energies up to 600 keV {x ~ 0.1). 
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IV. SUMMARY 



In summary, using a quasiclassical approximation we have obtained an expression for 
the dipole bremsstrahlung probabihty during a decay which is in agreement with the full 
quantum mechanical treatment of Papenbrock and Bertsch js] in a substantially wider region 
of the variable x = E^/Qa in comparison with all previous quasiclassical results. Our 
amplitude given in Eq. ([8]) contains both real and imaginary parts and, thus, includes the 
contribution from bremsstrahlung during the tunneling process. Our results demonstrate 
that the latter contribution is not negligible even for rather small x. As an illustration of 
these statements, we have considered the experimentally important case of ^^''Po. 

Furthermore, we find the quasiclassical expression for the contribution of the interference 
term of the dipole and quadrupole components to the double differential bremsstrahlung 
probability (with respect to the energy and the solid angle of the photon). This contri- 
bution turns out to be significant. Because of obvious limitations to the solid angle that 
can be covered by detectors in a realistic experiment, the angular distribution needs to be 
considered in the analysis of the experimental data, even though the quadrupole term makes 
a negligible contribution to the bremsstrahlung probability after integration over the entire 
solid angle. Using the expression for the dipole and quadrupole amplitudes presented in this 



work, the data analysis in our recent experiment was performed as described in Ref. |16| . 
and good overall agreement of our theoretical and experimental results was obtained (see 

n 

Fig. 5 of Ref. PJ[). We note, however, that a certain subtle question remains with respect to 
a next-to-leading order, nuclear mo del- dependent correction to the dipole-quadrupole inter- 
ference term, as discussed in Appendix |Bl These questions leave room for further interesting 
investigations in the context of under-the-barrier emission of bremsstrahlung in a decay in 
the future. 
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APPENDIX A: EFFECTIVE CHARGES 



The purpose of this appendix (see also 20|) is to clarify how the effective charges in 



Eqs. ([2]) and f[T^ for the dipole and the quadrupole terms arise in the interaction of a two- 
body system (charges eZ\ and 6^2, masses m\ and m2, e is the proton charge) with a photon 
of polarization e and wave vector q, as given by the part of the Hamiltonian corresponding 
to emission of a photon, 

Ri = e ■ pi exp (-1 q ■ Vi) 



rrii 
eZ2 



P2 exp {-iq- ra) . (Al) 



1712 

We define the total mass M = mi + m2, the reduced mass fi = mim2/M , the center-of- 
mass coordinate R = (mi/M)ri + {1712/ M)r2, and the relative coordinate r = Vi — r2. Let 
p and P be the momenta corresponding to the coordinates r and R, respectively. Then 
Pi = p+ {mi/M)P and p2 = —p+ {m2/ M)P . Writing the Hamiltonian (lAll) in the center- 
of-mass frame [P = 0) and performing its expansion over |q ■ ^ 1 up to the first term, 
we obtain 



Ht = —e 



z!S-iZiSq-r e-"^-", (A2) 



where the overall phase factor e""^ "'^ can be safely ignored. 
The effective charges are 



(1) _ Zi 1712 -Z2mi . . . 



(mi + m2) 



APPENDIX B: MATRIX ELEMENTS 



In this appendix we present some details of the derivation of our quasiclassical dipole ([H]) 
and quadrupole (fT3l) matrix elements. The wave function of the final state has the form 
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(see, e. g., §136 and §137 of [19|): 

— 5^i'e-^' {21 + 1) Rk,,i{r) Pi{n ■ A) , (Bl) 



2kf 

where n = r/r, A = kj/kj, Pi are the Legendre polynomials, and 

Rk,,i{^) = ^Fi{r^f,kfr) (B2) 
is the regular radial solution of the Schrodinger equation in a Coulomb field. 

ri 

When comparing to Eq. (5) of Ref. [8|, it is evident that the ansatz (lB2p for the final- 
state wave function corresponds to the neglect of the contribution from the irregular solution 
in the Couloumb field, which given by the term 2r~^Gi{'r]f, kfv) tana in the integral in 
Eq. (5) of Ref. ^. The basis for our approximation (]B2p is as follows. The asymptotics 
of the wave function at /cr ^ 1 is 2r~^ sin(/cr + fc + ^n), where fc corresponds to the 
asymptotics in a pure Coulomb field, and the phase shift ^at is due to the nuclear potential, 
which has a size of the order of tq. So, the asymptotic form of Rkj,i{r) is proportional to 
cos(57v)F/(?7/, kjr) + sm{SN)Gi{r]f, kfv). 

It follows from the quasiclassical approximation that the wave function at r ^ 
ref = Tt, where Vt = 2'r]f/{fiv) is the classical turning point, is given by Fi{rif,kfr) oc 



exp[—2^y2fivrif {^/rl — \/r)]. Therefore, 5^ oc exp[— 8?7/] ^ 1, and one is therefore led 
to the tentative conclusion that the term with the G function should be entirely neg- 
ligible. However, one might still object that the function G is exponentially large at 
r -C Tf, namely Gi{r]f,kfr) oc exp[2y/2/it^ — V^)]- order to convince ourselves 
that the contribution from the G functions is indeed small, let us consider as an ex- 
ample the term sm{6]\f)GQ{rii,kir) Gi{rif,kfr). This term is proportional to exp[— 4(?7j — 



rji)] exp(—Ay/2fIvf]^/r). Therefore, the main contribution to the transition matrix element 
from this term is given by the region r ~ l/{32firiv) = l/{32fi{zZa)) ^ tq. This contribu- 
tion is indeed small, and we can safely use the approximation ( 1B2|) for the final state of the 
a particle. 

The wave function of the initial state is given by an S state which consists of an outgoing 

wave at r ^ oo, 

z^,(r) = = -L- [Go{v^, Kr) +iFoim, hr)] . (B3) 

v47r v47rr 

14 



Substituting the wave functions into the transition matrix element and taking the integrals 
over the angular parts of the a particle wave functions, we obtain 



{f\Hj\t) 



2kf jj?uj 



e A 



+ Z^^e''' cose I dri?fc^,,2(r) 



ruj + 



Mr] 



(B4) 



where 6 is the angle between the vectors q and X, uj = \q\ = E^. Defining Ai and A/" as 



= -ir/i^/^^^ drRkf,i{r)Ro{r) 



2'Kk 



f Jo 



/ir 



(B5) 
(B6) 



and taking into account that the sum over the photon polarizations gives |e ■ Ap = sin^ 9, 
we obtain the double-differential bremsstrahlung probability as 



sin' ^ Id 



C =Z^^e'^'M + Z^^e'^Wcos^ ^ 



(B7) 



in agreement with Eq. fllOp . Then we use the quasiclassical approximation for the radial 
part of the wave functions (see, e.g., § 48, 49 of Ref. [l^). The matrix elements with the 
quasiclassical radial wave functions have been calculated using methods described in detail 



in Ref. 



211 ]. Although these methods are in principle well known, we present here some 



details of the calculation. 

Let us consider first the contribution Iq of the classically allowed region to the matrix 
element 

CO 

(B8) 



h = j drRkj,i{r)R^{r) 
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Using the standard quasiclassical wave function and assuming x <^ 1, we obtain 

oo 

J ^ L / ^''^ ^ig(sinh^9+i?) 

° Aif]J cosh2(^?/2) 



1 +^^l__^tanh(^9/2) 



1 

2f] 



i5(sinhi?+i9) 





oo 

i 



7] 



(I'd sinh.{-&) e 



^^e^^^^^"'^ , (B9) 







where bi = {I + 1/2)^ and the second expression is obtained after an integration by parts. 
According to 9|, the classically forbidden part can be incorporated by shifting the integration 
contour for ^ into the complex plane, via the replacement 'd ^ 'd + iir. As a result we arrive 
at 



(oo 
J dd sinh(79) e 



i5(smhi?-)9) 





oo 

I 



+ - j rf^9e'«(^*°'^''-''M . (BIO) 

Similarly, the leading in 1/f/ contribution to the quadrupole amplitude is determined by the 
integral [see the second term in square brackets in Eq. ( 1B6I) ] 



oo 

hi = j drRk^,2{r)ru;Ro{r) 



oo 



2f] 







The results ( 1B10|) and ( IBlip immediately verify Eqs. ([H]) and ( |T3l) . 

For the second contribution to the quadrupole amplitude, see Eq. (1B6I) . 

oo 

I22 = [ dri?fc^,2(r)— i?o(r), (B12) 


we find that this term is suppressed by a factor However, while J21 vanishes in the 

limit of a small photon energy x ^ 0, the contribution I22 tends to a nonzero constant 
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at a; = 0, i.e., even though I22 is parametrically suppressed by a factor 1/f/, it constitutes 
the dominant contribution to the dipole-quadrupole interference term at very small photon 
energies, due to its distinctive asymptotic behaviour. 

A precise calculation of the I22 contribution to the interference term is unfortunately 
hampered by the fact that the region of integration around r ~ ro gives an important con- 
tribution to the value of I22 due to the inverse power of r in the integrand, which thus 
makes the value of I22 nuclear model-dependent. In the data analysis of our recent experi- 
ment [1^, we therefore did not include the parametrically suppressed term l22- While good 
overall agreement between experiment and theory was obtained in this way (see Fig. 5 of 
Ref. [3]), the agreement is better at photon energies > 250 keV as compared to photon 
energies below this region. As the region of small coincides with the region where the I22 
term might contribute to the dipole-quadrupole interference term, its significance cannot be 
completely ruled out at present. 
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